Extinction of coherent backscattering by a disordered photonic crystal with a Dirac 

spectrum 
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Photonic crystals with a two-dimensional triangular lattice have a conical singularity in the spec- 
trum. Close to this so-called Dirac point, Maxwell's equations reduce to the Dirac equation for an 
ultrarelativistic spin-1/2 particle. Here we show that the half-integer spin and the associated Berry 
phase remain observable in the presence of disorder in the crystal. While constructive interference 
of a scalar (spin-zero) wave produces a coherent backscattering peak, consisting of a doubling of the 
disorder-averaged reflected photon flux, the destructive interference caused by the Berry phase sup- 
presses the reflected intensity at an angle which is related to the angle of incidence by time-reversal 
symmetry. We demonstrate this extinction of coherent backscattering by a numerical solution of 
Maxwell's equations and compare with analytical predictions from the Dirac equation. 
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Coherent backscattering is a rare example of an opti- 
cal interference effect that is systematically constructive 
in a random medium [l| . A reciprocal pair of waves (re- 
lated by time reversal symmetry) arrive in phase at the 
observer regardless of the path length or scattering se- 
quence. By measuring the angular profile of the diffu- 
sively reflected intensity for plane wave illumination and 
averaging over the random scattering, a peak is observed 
at a specific reflection angle 0, H, 3 . Under optimal con- 
ditions (same polarization of incident and reflected wave, 
transport mean free path I much larger than wave length 
A), the average peak and background intensity have ra- 
tio 7?. = 2 0, Q . This factor of two enhancement follows 
directly by comparing the coherent addition of intensi- 
ties (first sum the two wave amplitudes, then square to 
obtain the total intensity) with the incoherent addition 
(first square, then sum). 

Because coherent backscattering does not depend on 
the random scattering phase shifts, it is a sensitive tool 
to probe systematic phase shifts that contain information 
about internal degrees of freedom of the scatterers and 
the photons. One example in light scattering from cold 
atoms is the coupling of the photon polarization to the 
collective spin of the atomic ensemble, which can change 
the constructive into a destructive interference This 
effect is similar to the change from 7?. = 2 to 7?. = 1/2 in 
electronic systems with strong spin-orbit coupling [8| . 

A few years ago, Bliokh discussed an altogether dif- 
ferent and more dramatic switch from constructive to de- 
structive interference: ultrarelativistic fermions (with an 
energy much greater than the rest energy) would have a 
complete extinction of coherent backscattering, soTZ = 0, 
as a consequence of the Berry phase of tt accumulated 
along a closed trajectory by a half-integer spin pointing 
in the direction of motion. Indeed, it was previously no- 
ticed by Ando et al. in the context of graphene [l^l that 
the scattering amplitude s(0) for massless Dirac fermions 
vanishes when — > tt (with 4> the angle between the ini- 



tial and final wave vectors ki and kf). This absence of 
backscattering can be understood either in terms of the 
Berry phase difference of tt between time reversed scat- 
tering sequences or in terms of the antisymmetry 
S{ki kf) = —S{—kf —ki) of the scattering matrix 
of the Dirac equation Absence of backscattering 

in graphene might be measurable if scanning probe mi- 
croscopy can provide the required angular resolution of 
the electron flow [T3 |. 

Here we investigate an alternative realization of the ex- 
tinction of coherent backscattering that relies on photons 
rather than ultrarelativistic fermions. The half- integer 
spin required for the Berry phase of tt is produced by the 
Dirac-type band structure of a triangular-lattice photonic 
crystal [ll,[l3. By using photons rather than electrons 
the difficulty of angular resolved detection is avoided. We 
first discuss the effect at the level of the Dirac equation, 
and then test the predictions with a numerical solution 
of the full Maxwell's equations. 

We consider a photonic crystal with a two-dimensional 
(2D) triangular lattice structure. The hexagonal first 
Brillouin zone is shown in Fig. [1] Haldanc and Raghu 
[l3t showed that a pair of almost degenerate envelope 
Bloch waves (^"1, \E'2) = near a corner of the Brillouin 
zone can be represented by a pseudospin, coupled to the 
orbital motion. On length scales large compared to the 
lattice constant a and for frequencies near the degeneracy 
frequency ajD, the wave equation reduces to 

with Pauli matrices cr^,, ay. This is the 2D Dirac equation 
of a spin-i particle with zero mass and group velocity V£) 
(of order qlod)- The dispersion relation 

{u-u:D)^^vl{kl + kl) (2) 

has a double cone with a degeneracy at frequency lod 
(the so-called Dirac point). 
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FIG. 1: Hexagonal first Brillouin zone of a two-dimensional 
triangular lattice photonic crystal, with equifrequency con- 
tours centered at the corners. The three solid circles are re- 
lated by translation over a reciprocal lattice vector, so they 
are equivalent, while scattering from a solid to a dotted circle 
is suppressed if the scattering potential is smooth on the scale 
of the lattice constant a. The large dashed circle at the cen- 
ter is the equifrequency contour in air, included to indicate 
the refraction at the air-crystal interface (as expressed by Eq. 
©)• 



The Berry phase associated with the pseudospin degree 
of freedom was calculated in Ref. The solution of 

Eq. ^ with a definite wave vector k = {kx,ky) is 



{ui - ujd)/vd\ ^ 
kx + iky J 



cos(6i/2) 
e^f" sin(6'/2) 



(3) 



with C a normalization constant. The angles </>, 9 define 
the Bloch vector (cos </> sin 0, sine/) sin 6*, cos 6*), represent- 
ing the direction of the pseudospin on the Bloch sphere. 
Because of the dispersion relation ^ the angle 6 — 7r/2, 
so the Bloch vector lies in the x — y plane, pointing in 
the direction of k. The Berry phase (j>B is one half the 
solid angle subtended at the origin by the rotating Bloch 
vector. A rotation of fc by 360° in the x — y plane thus 
produces a Berry phase = n. 

Because of refraction at the interfaces a; = and x = L 
between the photonic crystal and air, we need to distin- 
guish the initial and final wave vectors ki, kf of enve- 
lope Bloch waves inside the crystal (velocity V]j) from 
the corresponding values fef"', kj" of plane waves outside 
(velocity c). For the crystallographic orientation shown 
in Fig. m the wave vectors before and after refraction (at 
a given frequency w) are related by [l^ 



+ 27r/3a, 



'■{k^"f 







(4a) 
(4b) 



with a the lattice constant and tj^) the frequency of the 
Dirac point. We denote k^" = Ifc^''! and k = |fc|. 

Before considering the diffuse reflection from the disor- 
dered photonic crystal, we address the specular reflection 
from the air-crystal interface that is present even without 
any disorder. A plane wave incident at an angle 

9 = arcsin(cA:f' '/a;) — arcsin(27rc/3tja + cki^y/cj) (5) 



FIG. 2: A pair of reciprocal waves that interfere destruc- 
tively, resulting in the extinction of coherent backscattering. 
Arrows at the air-crystal interface indicate the incident plane 
wave (red solid lines), the diffusively reflected wave (blue solid 
lines), and the specularly reflected wave (black dashed lines). 
The two initial waves (wave vector kf" , refracted to ki) follow 
time-reversed sequences of scattering events (dark circles) , to 
end up in two final waves (wave vector kf, refracted to kj") 
with a phase difference of TT -I- (fci -I- fe/) -(Jr. Forfe/ — — fci only 
the Berry phase difference of tt remains. As a consequence 
of this destructive interference, the intensity of the reflected 
wave indicated in blue is suppressed. By measuring the an- 
gular profile of the average reflected intensity a minimum of 
nearly zero intensity will result at this angle. 



is specularly reflected at an angle ^spoc — t^^9. The inter- 
face reflectivity i?int (the fraction of the incident photon 
flux that is specularly reflected) follows from the transfer 
matrix of the air-crystal interface calculated in Ref. [1^1 . 
In the approximation of maximal coupling we find 



-Rint(fcy) = 



1 - v/1 - {ky/kf 



1 + VI - {ky/kf 



(6) 



hence the interface reflectivity is zero for ky — 0. Numer- 
ical solutions of Maxwell's equations give i?int(0) ~ 
0.05 for Lu near cod, so this is a reasonably accurate ap- 
proximation. 

Disorder inside the photonic crystal produces a back- 
ground of diffusively reflected waves in an angular open- 
ing 59 ~ k/k^" ^ 1 around 6*0 = tt — arcsin(27rc/3wa). 
The reciprocity angle 0* for coherent backscattering is 
related to the incident angle 9 of Eq. ^ by 



TT — arcsin(27rc/3wa — cki^y/uj). 



(7) 



We will choose ki^y small (<C fc) but nonzero, so that 
^int «C 1 while still coherent backscattering at angle 9^, 
can be resolved from specular reflection at angle ^spoc. 

The extinction of coherent backscattering by destruc- 
tive interference of reciprocal waves is illustrated in Fig. 
[21 The two series of time reversed scattering events 
S+ = fcj — > fci — > fe2 — * • • • — > fe„_i — > kji ~~ki and 
>S'— = ki > k^ > kji—i > ■ • • > k2 > ki > ki 
have the same scattering amplitude up to a Berry phase 
difference 0_b = tt [1, . This destructive interference 
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suppresses the reflected intensity at angle 6^, . If the final 
wave vector k f deviates from the exact backscattering di- 
rection — fcj, the phase difference = (j)B + {ki + kj)-Sr 
depends on the separation 5r of the first and last scat- 
tering events The Berry phase difference (j)B remains 
equal to tt, because both the Berry phases accumulated 
along 5+ and 5*- are incremented by the same amount 
(half the angle between fc/ and —ki). 

By including the Berry phase in the theory P, Q 
of coherent backscattering of scalar waves in the weak 
scattering regime {l^ X), it follows that the incoherent 
part i?o of the reflectivity R remains unaffected while the 
interference part 5R acquires an overall factor cos 0b Q: 



R = Rq + SR cos (f)B- 



(8) 



The entire angular profile of coherent backscattering in 
the Dirac equation (where (/)b = tt) may therefore be 
obtained, for / 3> A, by simply changing the sign of the 
known results for SR for scalar waves (where 0s ~ 0). 
In particular, for ki^y ^ k and Sk = kf^y + ki^y <C 1// one 
has 



Rn 



1 



(1 + WO 1 



L + 2zo 



(9) 



X {l - (/ + zo)(5fccoth[(5fc(L-|-2zo)]}. (10) 

The normalization factor N ~ kW/2TT is chosen such 
that R is the fraction of the incident photon flux which is 
reflected in a single transverse mode when kf^y = 2'im/W 
{n = 0, ±1, ±2, . . .) is discretized by periodic boundary 
conditions a.t y ~ Q and y = W . (This normalization 
is chosen to simplify the comparison with the numerical 
calculations described later on.) 

Eqs. Q and (fTO|) are approximate results from the 
radiative transfer equation, accurate for a disordered slab 
of thickness L not much smaller than the transport mean 
free path The parameter zq, the so-called extrapolation 
length, depends on the reflectivity i?int of the interface 
between the photonic crystal and vacuum, according to 

M 

zq = \ttI- C*! = 7 / Rint{ky)dky, (11a) 

i — C-i Jo 

C2 = - [ jl-{ky/kYR,^,{ky)dky. (lib) 

Substitution of Eq. ^ gives the simple answer zq ~ I for 
the extrapolation length of the air-crystal interface. 

Collecting results, we arrive at the following line shape 
of the reflectivity near the reciprocity angle: 



4 / I 

R= lSkcoth\dk(L + 21)] 

nN \ ^ ^ L + 21 



The suppression p2|) of the reflected intensity within a 
narrow angular opening 50 ~ 1/kl around the reciprocity 
angle 0* is compensated by an excess intensity 6R ~ 1/fc/ 
of the diffusively reflected wave at angles away from 6**. 
This compensation is required by current conservation 
[l3|, but difficult to observe for kl ^ 1. We will therefore 
not consider it in what follows. 

We now compare the analytical predictions from the 
Dirac equation with a numerical solution of Maxwell's 
equations. As in earlier work [1^1, we use the meep soft- 
ware package llSl to solve Maxwell's equations in the 
time domain [l9[ for a continuous plane wave source 
(time dependence cx e^"*) switched on gradually. We cal- 
culate the reflected intensity, projected onto transverse 
modes, as a function of time and take the large time limit 
to obtain the stationary reflectivity R. We took the time 
sufficiently large that the sum of the total transmission 
and total reflection differs from unity by less than 0.03. 

The photonic crystal consists of a triangular lattice 
of parallel dielectric rods (dielectric constant 8.9, radius 
Pq = 0.3 a) in air. The orientation of the lattice is as 
shown in Fig. ^ The magnetic field is taken parallel 
to the rods (TE polarization). The conical singularity 
in the band structure is at frequency lod = 3.03 c/a, 
with a slope dui/dk = vb ^ 0.432 c. We discretize 
the transverse wave vector by means of periodic bound- 
ary conditions: ky = 2'!rn/W, n = 0,±1,±2,..., with 
W = 1200 a. The longitudinal dimension of the lattice 
is taken at L = 60^/3 a (corresponding to 121 rows of 
dielectric rods). The frequency lo of the incident plane 
wave is chosen such that lu — ujd = 0.206 c/a is sufficiently 
large that N = 88 ^ 1, but sufficiently small that the 
trigonal distortion of the circular equifrequency contours 
is insignificant. 

Disorder is introduced by randomly varying the radius 
p{r) of the dielectric rods at position r, according to 
p{r) = Pa + Sp{r). The variation Sp{r) is spatially cor- 
related over a length ^ in a Gaussian manner. 



(13) 



3nN 



{LI + 2f)Sk^ for 6k 0. 



(12) 



where {i = 1,2,.. .A/") is a randomly chosen point in 
the crystal and Ai is a randomly chosen amplitude in 
the interval (-A, A). We took TV = 3200, A = 0.04 a, 
^ = 2a. It is essential that the correlation length ^ of the 
disorder is larger than the lattice constant a to minimize 
scattering between the two inequivalent corners of the 
Brillouin zone (solid and dashed circles in Fig. [1]). For 
these disorder parameters we found that only about 2% 
of the incident flux is reflected into the opposite corner. 

The numerical data in Fig. [3] is for the incident wave 
vector ki^y = —32tt/W, averaged over 80 realizations of 
the disorder. The angle of incidence (measured relative 
to the positive x-axis) is 9 = 0.67 rad = 38.4°. The 
specularly reflected wave at an angle ^spcc ~ — 9 ~ 
2.47 rad = 141.5° has transverse mode number n = —16. 
The reciprocity angle ^ is 9* = 2.40 rad = 137.5°, cor- 
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FIG. 3: Plot of the reflectivity R (fraction of incident flux 
reflected in a single transverse mode) as a function of the 
transverse mode index n (related to the transverse wave vec- 
tor by kf^y = 2nn/W). The specularly reflected flux is at 
n = — 16 and the extinction of coherent backscattering is at 
n = +16 (indicated by arrows). The curve is calculated nu- 
merically from Maxwell's equations. The dashed and dotted 
lines are, respectively, the analytical predictions (|9]) and (|12p 
for the incoherent background and the line shape near the 
extinction angle, for a transport mean free path of Z = 10 a. 



responding to the mode number n = +16. 

Both the specularly reflected wave at n = — 16 and the 
extmction of the coherent backscattering at n = +16 are 
clearly visible in Fig. [31 The extinction at the reciprocity 
angle is not complete, presumably because of scattering 
between inequivalent corners of the Brillouin zone. The 
analytical theory predicts a parabolic line shape of the re- 
flectivity around the reciprocity angle, given by Eq. ([7]), 
the width of which depends on the value of the trans- 
port mean free path /. A fit to the numerical data gives 
I = lOa (dotted curve in Fig. [3]), which then implies 
a background of incoherent diffuse reflection at a value 
Rq ~ 0.006 which is somewhat larger than the numerical 
data (dashed horizontal line). 

In conclusion, we have shown that coherent backscat- 
tering of radiation from a disordered triangular lattice 
photonic crystal is extinguished at an angle that is re- 
ciprocal to the angle of incidence. This effect has an 
analogue for electrons d, [l3| , but it should be easier to 
observe for photons because angular resolved detection 
is more feasible. The observation of the extinction of co- 
herent backscattering would be a striking demonstration 
of a spin- 1/2 Berry phase in a disordered optical system. 
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